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Professors  Birnbaum  and  Neyman  have  elected  me*  (as  they 
assert  my  native  language  Is  English,  an  honor  rarely  accorded  to 
Australians),  to  thank  the  organizers  and  to  suggest  that  this 
dialogue  between  biologists  and  statisticians  be  continued  In 
future  meetings.  It  Is  painfully  clear  that  we  have  only 
scratched  the  surface  of  the  problem.  This  problem  Is:  how  to 
understand  environmental  hazards  to  health  through  experimental 
and  epidemiological  research,  and  how  to  provide  wise  advice  to 
those  who  roust  frame  policies  for  regulation.  It  Is  truly  an 
enormous  task. 

Oak  Ridge  Is  a natural  place  for  such  a meeting  since  much  of 
the  environmental  problem  arises  from  the  demand  for  more  energy 
and,  of  course.  Oak  Ridge  owes  its  existence  to  nuclear  energy. 

The  most  studied  hazard  Is  radiation  and  the  largest  animal  ex- 
periments have  been  conducted  here.  But  chemicals  are  an  Increas- 
ing menace  to  health.  In  particular,  the  use  of  fossil  fuels 
to  generate  energy  raises  many  medical  problems.  Thus  the 
Institute  for  Energy  Analysis  at  Oak  Ridge  Is  a natural  host 
for  this  meeting.  Of  course,  such  work  Is  being  done  by  many 


groups,  e.g.,  the  National  Research  Council  has  a committee  con-  j 

sidering  the  "Research  Needs  for  the  Health  Effect  of  Fossil  Fuels."  1 
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1.  INTRODUCTION  - COMMENTS  ON  THIS  CONFERENCE 
I have  been  asked  to  make  my  observations  on  this  meeting. 

The  papers  have  reviewed  some  literature*  discussed  some  experiments 
and  their  analysis*  some  Ideas  have  been  suggested  that  are  worth 
following  up*  and  many  essential  things  have  not  been  mentioned. 

Speakers  have  continually  referred  to  cancer.  This  is  no^  the 
greatest  public  health  problem!  I believe  that  the  end  of  live  Is 
less  Important  than  the  quality  of  life*  which  means  that  I think 
\ chronic  diseases  should  have  been  the  main  topic.  When  wc  leave 

I radiation  and  turn  to  chemicals  (e.g.*  the  effects  of  effluents  from 
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1 a fossil  fuel  power  station)*  we  think  of  debilitating  thincjs  such  as 

asthma*  etc.  And  yet  everyone  has  been  concerned  with  the  causes  of 
death  and  the  length  of  life. 

r Below  I have  given  my  views  on  a number  of  Issues.  Ky  major 

I points  are  - 

(1)  there  are  n£  competing  risks  In  the  ILLNESS  - DEATH  MODEL 
(11)  it  seems  very  hard  to  frame  a biologically  reasonable 
model  to  make  sense  of  the  potential  lifetime  concept* 
but  an  attempt  Is  made  In  I6. 

2.  CAUSES  OF  DEATH 

For  human  populations  the  causes  of  death  usually  are  very 
hard  to  determine.  We  have  seen  that  even  with  the  most  lavish  | 

studies  of  mice*  they  are  equally  vague.  Thus  Nr.  Neyman  has  very 
sensibly  tried  to  persuade  us  that  we  should  try  to  work  with  what 
Is  a little  better  determined — the  medical  description  of  the  state 
of  the  deceased  mouse.  For  humans  who  have  not  been  under  con- 
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tinuous  good  medical  care»  even  this  will  be  vague,  too.  With 
computers,  such  complex  Information  can  be  studied.  More  empirical 
data-analytl cal  methods  (see,  e.g.,  the  new  book  by  Gnanadesikan) 
will  have  to  be  used. 

The  statisticians  present  learned  about  the  great  contrasts  be- 
tween animals  (usually  mice)  and  humans — the  permissible  errdyet  of 
exporlments  and  the  causes  of  death.  The  extrapolation  from  a.ii1ma1 
experiments  to  man  Is  agreed  to  need  a real  understanding  of  ::zuses 
and  mechanisms. 

3.  COMPETING  RISKS 

Much  of  our  time  was  spent  on  "The  Competing  Risks"  prcLlem  in 
different  ways  and  at  different  mathematical  levels  and  alway;  with 
different  notations!  The  basic  problem  seems  clear  with  no  .7.2ihema- 
tlcs  at  all.  Let  us  assume  that  there  Is  a list  of  C (primary) 
causes  cf  death.  Classify  and  count  all  deaths  In  a large  human 
population  In  a year  and  divide  each  number  by  the  population  size. 

If  the  first  cause  of  death  Is  eliminated.  It  seems  obvious  that, 
with  no  further  assumptions,  there  Is  no  way  of  predicting  the  C-1 
ratios  that  will  be  seen  In  the  future.  We  simply  don't  knew  what 
will  be  the  fate  of  those  who  formerly  died  from  the  first  cause. 

Biologically,  we  know  that  the  susceptibility  of  an  Individual 


to  various  pathological  states  (as  well  as  eye  colors,  etc.)  will 
depend  upon  his/her  genotype.  The  Incidence  of  these  conditions 
will  depend  upon  the  life  history  of  the  Individual,  I.e., 


the  Individual's  environment  In  the  largest  sense  of  the  word.  In  a 
population  there  Is  variation  In  both  genotype  and  environment.  At 
least  In  principle  these  two  sources  of  variation  would  have  to  deter- 
mine the  Joint  distribution  of  the  potential  times  to  death  the 

various  causes  of  an  individual  chosen  at  random  from  this  population. 
If  we  can  make  sensible  definitions  of  these  potential  lifetimes 
(see  56).  Now,  the  genetic  component  is  weaker  In  some  conditions 
than  In  others.  If  also  there  Is  1 1 ttle  correlation  betwcsi:  Indivi- 
dual environments,  we  may  expect  "potential  life  times"  for  these  con- 
ditions to  be  roughly  Independent.  If  the  genetic  background  Is  homo- 
geneous (e.g.,  an  Inbred  mice  strain),  we  need  only  worry  sbeut  the 
environmental  correlations.  Conversely,  If  the  environment  1c  con- 
stant, only  the  genetic  variation  need  be  considered. 

Cve»‘yonc  knows  of  the  difficulties  In  the  Nature-Nurture  argu- 
ment---!.Q.  and  race  Is  a notorious  example.  Even  Identical  twin 
studies  may  be  criticized.  The  genetics  of  quantitative  characters 
Is  simply  not  in  a state  to  help  us. 

The  several  paragraphs  above  were  written  In  the  hope  that  they 
will  throw  some  biological  light  on  the  statistician's  discussion 
of  "Independence"  and  "dependence"  In  relation  to  "Competing 
Risks."  Recall  also  that  the  notion  of  a cause  of  death  Is  easily 
criticized.  Thus  the  statisticians'  formulation  of  the  problems, 
however  amusing  It  may  seem  to  them,  may  not  be  very  good.  That 
It  leads  to  conundrums  like  "non-ldentiflablllty"  may  not  be  a 
rsi:s&  foY*  tears.  The  problem  Is  real  enough  however.  And  yat  i 


have  never  been  able  to  convince  medical  people  that  It  might 
upset  their  conclusions  that  a treatment  was  reducing  a death  rate» 
for  example!  I have  the  feeling  that  this  problem  may  be  over- 
emphasized.  The  theory  appears  also  under  the  subject  heading  of 
Reliability  Theory.  However,  I am  convinced  this  latter  theory 
has  not  made  any  technical  devices  more  reliable,  though  many  of 
us  have  enjoyed  writing  and  reading  papers  about  It. 

I win  return  to  this  argument  In  i6. 

4.  ON  HERETICS  AND  HEROES 

I am  personally  attracted  to  people  who  rationally  advocate 
non-estabi Ishmsnt  views.  Two  come  to  mind  In  this  general  area 
and  I mention  them  In  case  they  are  unfamiliar  to  some.  P.  R.  J.  Burch 
has  written  two  books  and  many  papers  In  which  the  ages  of  Incidence 
or  onset  of  many  conditions  has  been  studied.  He  seems  to  have  done 
more  work  In  this  area  than  anyone  else.  In  a recent  letter,  he 
charmingly  described  himself  as  a "recognized  heretic."  While  well 
aware  of  the  many  difficulties  with  medical  data,  he  Is  not  unduly 
worried  about  his  curves  being  upset  by  competing  risks. 

I turn  now  to  cancer  to  make  several  points.  In  this  area, 

Arrain  Braun  Is  something  of  a heretic.  His  book.  The  Biology  of 
Cancer,  makes  very  good  reading.  He  does  not  accept  the  almost 
universal  equating  of  mutation  and  cancer.  (See  also  my  note  In 
P.N.A.S.,  April  1977).  Certainly  all  speakers  at  this  meeting  made 
this  9«5umpt1on---of  course  It  may  well  be  true  In  these  specific 

L'._ 


cases--osteosarcoma  due  to  bone-incorporated  strontium  90  and 
melanoma  due  to  U.V.  In  sunlight.  Other  cases  where  It  Is  an 
attractive  hypothesis  are  the  early  childhood  cancers  like  re- 
tinoblastoma. The  models  presented  have  a natural  appeal  to  an 
applied  mathematician/statistician.  But  they  seem  to  raise  as 
many  questions  as  they  answer.  (Of  course,  this  is  characteristic 
of  any  good  theory.)  One  point  that  worries  me  Is  the  following. 

In  the  resting  state  the  DNA  In  a nucleus  Is  a supercollsd, 
tangled  mass  like  a full  pot  of  spaghetti.  To  take  the  Groer-Marshal 1 
theory,  an  a particle  going  through  this  mess  Is  likely  to  do 
a lot  of  damage.  Some  of  this  damage  Is  supposed  to  lead  to  a 
mutation  which  Is  a step  to  carcinogenesis,  some  Is  to  be  repaired 
and  some  to  kill  the  cell.  It  seems  Intuitively  to  me  that  some 
damage  will  be  Irreparable,  and  some  will  be  repaired,  often  with  error 
so  that  many  mutations  will  be  produced.  The  progeny  of  some 
mutated  cells  will  not  survive  the  competition  with  normal  and 
other  mutated  cells.  Hany  mutations  may  be  synonomous  or  selectively 
neutral  and  so  for  growth  purposes  be  equivalent  to  normal  cells. 
Finally  some  mutated  cells  are  freed  from  growth  restraints  and 
become  cancer  cells.  Has  anyone  looked  for  mutants  among  sur- 
viving cells  and  sought  to  relate  their  prevalence  to  dosage? 

For  they  should  obey  the  same  laws  as  cancer-transformed  cells.  Such 
mutants  should  be  detectable  electrophoretical ly. 

There  Is  a large  literature  no  doubt  examining  critically 
the  relative  sizes  of  DNA  molecules  and  their  radiation  cross- 
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sectlons.  It  Is  a little  dangerous  to  think  about  cloud  chamber 
tracks  which  are  visible  to  the  naked  eye  and  make  deductions  about 
molecules . 

Since  It  Is  useful  to  read  Into  the  record  of  conferences 
references  to  highly  relevant  work,  I would  like  to  mention  the 
P.N.A.S.  1976  paper  by  Neyman  and  Purl  on  a model  for  radiation 
effects  and  a 1974  book  Intrinsic  Mutagenesis  by  Sir  Macfarlane  Burnet 
which  Is  concerned  with  the  whole  problem  of  this  conference. 

5.  ILLNESS  AHD  DEATH  MODELS 

These  models  stemming  ultimately  again  from  Hr.  Neyman,  have 
been  discussed  by  C.L.  Chlang  here  and  In  his  well-known  book. 

In  passing  the  recent  result  of  P.  Clifford  (P.N.A.S.,  1977)  was  al- 
luded to.  This  Is  such  an  Important  comment  on  such  models  that  It 
Is  a pity  It  was  not  discussed  at  length.  Hence  I will  give  a brief 
summary. 

The  random  progression  of  the  health  state  of  an  organism  might 
well  be  a Markov  process  If  the  description  of  the  state  Is  com- 
plex enough.  (There  was  some  discussion  of  this  and  It  seems  that 
In  practice  the  state  designation  Is  too  simple  for  this  to  be 
valid.)  There  will  be  many  transition  functions  to  be  estimated 
In  any  actual  experiment.  In  assessing  the  effect  of  a changing 
environment  we  will  be  particularly  Interested  In  transition  func- 
tions that  are  time  dependent. 


If  one  considers  the  simplest  conceivable  such  model,  Clifford 


shows  that  the  largest  of  survival  and  sacrifice  experiments  will 
lead  to  ambiguity!  Let  H stand  for  healthy,  for  having 

disease  1,(1»1.2)  and  *12  for  having  both  diseases.  Suppose 

only  forward  transitions  are  possible  and  that  death  can  occur  In 
any  of  the  four  live  states.  The  following  figure  should  now  be 


If  (f)  probabilities  of  being  In  live 

n 5 j 

states  H,  D,  etc.,  at  time  t,  then  with  a dot  for  d/dt. 


r 


-9- 

Denoting  the  death  states  by  0^^,  Dj,  0^,  0^2*  similarly 


(t) 


12 


<2) 


Now  with  an  enormous  (i.e.,  "gedanken")  experiment  and  observ 
ing  and  sacrificing  at  very  fine  time  intervals,  we  can  determine 
all  the  P's  as  functions  of  time  so  that  aquations  (2)  determine 
the  y's.  Writing  + VjjP|,  » Rj,.  etc.,  which  are  now  known 


functions  of  time,  (1)  becomes 


“ •'h " 

1 
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1 
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These  equations  may  be  written  R *>  PV  where  P clearly  has  a 
zero  determinant.  Thus  there  is  no  unique  solution  for  the  y's, 
although  Vj  + V2  is  determined. 

The  "potential  life  time"  model  asserts  that  with  Causes  of 
death,  there  is  a joint  distribution  of  potential  times  to  d.ith 
^1'  death  T is  given  by 


T » min(Tj T^.)  . (4) 

This  makes  no  sense  whatever  in  this  model.  But  even  more  striking, 
there  are  no  competing  risks  in  this  model!  For  when  the  individual 
is  in  H,  a "healthy  death"  is  the  only  cause,  when  in  Sj,  this  is 
the  only  cause  of  death,  etc.  Thus  if  we  follow  Mr.  Neyman's  sugges- 
tion of  giving  a pathological  description,  and  there  are  no  back 
flows,  these  models  shed  no  light  whatever  on  competing  risks. 

6 . COMPETING  RISKS  AGAIN 

I would  like  to  develop  the  genetical  picture  in  53  to  give  yet 
another  account  of  the  competing  risk  problem.  My  excuse  is  that  I 
believe  it  throws  some  new  light  on  two  aspects  of  the  issue,  (1)  the 
definition  of  "potential  life  times,"  (ii)  the  question  of  dependence. 
This  work  was  done  after  the  conference  and  for  most  of  this  period 
I v/as  convinced  that  current  formulation  involving  (4)  made  no  sense 
whatever! 

For  any  individual  i,  we  suppose  that  the  chance  that  he  dies 

from  cause  c in  (t,t  + dt)  is  y(t,c,i)dt  where  c » 1,2 C. 

If  only  cause  C were  operating,  he  would  have  a lifetime  T^(1). 
say,  where 

(1)  Prob  (T^(i)e(t,t  + dt))  = exp^-/pP(t',c,i  )dt'|p(t,c,i)dt 

(2)  Prob(T^(i)  > t)  » exp|-/^y(t' ,c,i )dt' j . 

We  call  T^(i)  his  potential  lifetime  for  cause  of  death  c.  They 

are  clearly  non-observable  since  several  causes  act  simultaneously. 

When  all  C causes  are  competing  for  his  life,  define 

C 

=*  z M(t,c,i)  . 
c»l 


(3) 


y(t,i ) 


Then 


(4) 


P(  1 dies  of  cause  c at  an  age  in  (t,t  + dt)) 

= exp|-/Ju(t' ,i)dt' ju(t,c,i)dt  , 

C 

(5)  = n Prob(T  .(i)  > t+dt)Prob(T  (i)e(t,t+dt))  , 

c»l 

c'fc 

(6)  - Prob[T^(i)  = min(T.{i)....,T^.(i)J  ,T^(i)£(t,t+dt)| . 


(5)  follows  from  (4)  by  (1)  and  (2)  on  observing  that  we  need  only  the 
term  of  order  1.  (5)  and  (6)  make  it  clear  that  with  this  model  the 

individual's  potential  lifetimes  may  be  regarded  as  independent  and  he 
dies  at  the  age  of  the  least  one.  Further,  the  chance  of  death  for 
i in  {t,t+dt)  is 

exp|-/pu(t' ,i  )dt' jp(t,c,i  )dt  = exp|-/Q^p(t' ,i)dt' jy(t,i) . 

Setting  1 - F^{t)  * Prob(  i living  longer  than  t),  we  have 

(7)  1 - F^(t)  - exp[-/Q%(t' .i)dt'j  . 


Now  suppose  that  the  population  contains  a fraction  p^  of  in- 
dividuals like  i.  For  example,  this  group  all  have  the  same  geno- 
type and  environmental  history.  Then  the  life  distribution  in  the 
population  is 


(8) 


F(t)  = Ep^F^(t) 

<=  1 - Ep^  exp^-/p  w(t' ,i  )dt' j 
C 

= 1 - Ep.  n P(T.(i)  > t) 

’ c=l  ^ 


and  probability  density  of  life  spans  terminated  by  cause  c»  say, 

at  t is  the  average  of  (4)  or  (5)  or  (6), 


(9) 


P(L^c(t,t+dt))  = n(t,c,1  )exp|-/^^y(t' ,1)dt'jdt 


while 

(10) 


P(TJi)e(t.t+dt))  n P(T..(i)  > t)dt 
1 ^ c'fc 


P(Lc>  t)  * j"P(L^e(t".t"+dt“)). 


Thus  one  cannot  treat  the  L 's  as  Independent  variables  and  use 

c 


their  mlnijaum  as  we  could  the  T (i)'s.  The  L *s  could  conceivably 

C w 

be  called  the  potential  lifetimes  of  a randomly  chosen  individual.  If 
t-hs  "acting  alone  aspect"  is  valued  we  might  consider  instead 
where 

Prob(D^c(t,t+dt))  = Ep^  w(t,c,i)exp|-J^M(t' ,c,i)dt'J 

Prob{0^>t)  =»  Ep^  exp^-|pp(t' ,c,i  )dt'|  . 

These  do  not  seem  to  be  good  candidates. 

If  one  cannot  classify  individuals,  a detailed  study  of  the  popu- 
lation will  only  reveal 
-1 


(11)  f(dt)"^*fraction  of  individuals  ] ^ n » u(t  c )sav 

[ dying  of  c in  (t.t+dt)J  ^^i  ’^vt.c.ij  wit,c,i,say. 


(12)  f(dt)"^*fraction  of  individuals 
dying  in  (t,t+dt) 


I » E w(t,c),  * w(t),  say. 
^ c 


From  (12)  we  will  assert  that  the  life  distribution  is 

(13)  F(t)  = 1 - exp[-/^V(f)dfj  . 

There  is  always  a function  ir( • ) such  that 

(14)  exp^-j  ir(t)dtj  » Ep^  exp|-|  y(t',i)dt'| 


n 

S I 


j 
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as  may  be  seen  explicitly  by  taking  logarithms  and  differentiating. 
Similarly,  we  will  write,  instead  of  (9), 


exp  ir(t'  ,c)dt' ju(t.c)  . 


It  is  clear  that  even  if  the  p^  are  known  (but  not  which  individuals 
belong  to  type  i),  knowledge  of  ir(t,c)  does  not  determine  the 
y(t,c,i). 

With  this  set  up,  suppose  cause  C is  eliminated.  This  means 
that  p(t,C,i)  = 0,  for  all  i,t.  Equations  (1)  to  (7)  change  very 
simply.  For  example,  (7)  now  reads 

c-1 

(15)  1 - Ff(t)  « n P(T^(i)  > t)  , 

’ c=l  ® 

so  that  (8)  now  reads 


(16) 


1 - F*(t)  - zp.  n P(T  (i)  > t)  . 
c»l 


Thus  (16)  cannot  be  determined  unless  the  individuals  can  be  classified. 
If  they  can  be  classified,  all  the  problems  disappear  with  a 
sufficiently  large  experiment. 
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